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Abstract 

We present the first study of treating b, c, and s quarks as Dirac 
fermions in lattice QCD with exact chiral symmetry. For 100 gauge 
configurations generated with single-plaquette action at (3 = 7.2 on 
the 32 3 x 60 lattice, we compute point-to-point quark propagators for 
33 quark masses in the range 0.01 < m q a < 0.85, and measure the 
time-correlation function of (pseudo-)scalar, (pseudo-)vector, and ten- 
sor mesons, for symmetric and asymmetric quark masses respectively. 
The lowest-lying mass spectra of mesons with quark contents bb, cb, 
sb, and cc are determined, together with the pseudoscalar decay con- 
stants. Our results are sumarized in Tables 1-5. Furthermore, we 

also determine the b and c quark masses, mj^^(m(,) = 4.65(5) GeV, 

m^ S (m c ) = 1.16(4) GeV. 
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1 Introduction 



Currently, one of the most important themes in high energy physics is to find 
out whether there is any new physics beyond the Standard Model (SM). In or- 
der to identify any discrepancies between the high energy experimental results 
and theoretical values derived from the SM, the latter have to be obtained in 
a framework which preserves all crucial features of the SM. Otherwise, it is 
difficult to determine whether such a discrepancy is due to new physics, or just 
the approximations (or models) one has used. 

So far, the largest uncertainties in the theoretical predictions of the SM 
stem from the sector of the strong interaction, namely, QCD. Theoretically, 
lattice QCD is the most viable framework to tackle QCD nonperturbatively 
from the first principles. However, in practice, any lattice QCD calculation 
suffers from the discretization and finite volume errors (which can be system- 
atically improved). 

Since all quarks in QCD are excitations of Dirac fermion fields, it is vital 
to preserve this important feature in any approaches to unveil new physics 
beyond the SM. Theoretically, the most appealing lattice fermion scheme is 
the Domain- Wall/Overlap fermion [1, 2, 3], which preserves the exact chiral 
symmetry at finite lattice spacing [4], thus provides a proper formulation of 
QCD on the lattice. However, in practice, it is difficult to accommodate all 
quark flavors with presently accessible lattice sizes, since the quark masses 
span at least 3 orders of magnitude (from m u ~ 5 MeV to mj ~ 4600 MeV), 
even excluding the top quark. In fact, one still cannot put the physical u/d 
quark on a lattice, with sufficiently large number of sites in each direction such 
that the discretization and finite volume errors are well under control. Thus, 
one can only perform simulations at unphysically heavy u/d quark masses, 
and then use chiral perturbation theory to extrapolate to the physical limit 
(i.e., m n ~ 140 MeV). This scenario is not expected to have dramatic changes 
in the next few years. 

Nevertheless, if one only aims at putting b, c, and s quarks on the lattice 
(with m q a < 1, and M^L > 4), then a lattice of size ~ 32 3 x 60 (with inverse 
lattice spacing a^ 1 = 7.68 GeV) seems to be sufficient for this purpose. As we 
will see below, the meson mass spectra turn out in good agreement with the 
experimental values, for quark contents bb, cb, cc, and sb. 

Note that our approach is fundamentally different from other lattice QCD 
calculations using relativistic heavy quark formalism, or heavy quark effective 
theory (HQET), or non-relativistic QCD (NRQCD), in which the systematic 
errors are difficult to control. 

In this paper, we will restrict ourselves to mesons with quark contents 
bb, cb, cc, and sb. Our results of the masses and decay constants of the 
pseudoscalar mesons B s and B c , and also the masses of the vector mesons B* 
and B*, have been presented in Ref. [5]. For the mass spectra of spin- 1/2 and 
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spin-3/2 baryons with quark contents bbb, bbc, bcc, bcs, bss, and ccc, our 
results will be presented in a forthcoming paper. 



2 Lattice quarks with exact chiral symmetry 

To implement exact chiral symmetry on the lattice [1, 2, 3, 4], we consider the 
optimal domain- wall fermion proposed by Chiu [6, 7]. The action of optimal 
domain-wall fermion can be written as [7] 

N s +1 

A F = Yj Y^{x,s){(uj s D w (x,x') + 5 X)X i)5 ssI 

s,s'=0 x,x' 

+(uj s D w (x,x') - 5 x ^)(P + 5 slyS - 1 + P_5 s / )S+ i)}^(x',s / ) 

with boundary conditions 

P + ip(x, — 1) = — r m q P + if)(x, N s + 1), 

P-ip(x, N s + 2) = — r m„ P-ip{x, 0), r = , 

2m 

where m q is the bare quark mass, and the weights {u s , s = 1,---,N S } are 
specified by the exact formula derived in Ref. [6] such that the system possesses 
the maximal chiral symmetry for any fixed N s . Here H w = 75-D™, and D w is the 
standard Wilson Dirac operator plus a negative parameter —mo (0 < mo < 2). 
The quark fields are constructed from the boundary modes at s = and 
s = N s + 1 with u = uipf a+ i = [7]: 

q(x) = ^[P^(x,0) + P + ij(x,N s + l)], 
q(x) = y/r Up(x,0)P+ + xp(x,N 8 + l)P. 

After introducing pseudofermions with m q = 2mo, the generating functional 
for n-point Green's function of the quark fields can be derived as [7], 

_ /[dC7le-^det[(D c + m q ){\ + rD c y 1 } exp { J(D C + m q )- l j) 
Z ^ J] ~ J[dU]e- A 9det[(D c + m q )(l + rD.)- 1 ] (1) 

where A g is the action of the gauge fields, J and J are the Grassman sources 
of q and q respectively, and 

1 - 75 3 opt 

q 1 - nf=! Ts 



1 - uj s K 



w 



1 + oj s H w 
2 



From (1), the valence quark propagator in background gauge field is 

5 2 Z[J, J] 



(q(x)q(y)) 



(D c + m g ) 1 



Ux,y 

J=J=0 



5J(x)5J(y) 

where D c is exactly chirally symmetric (D c , y 5 +'-f 5 D c = 0) in the limit N s — > oo, 
and its deviation from exact chiral symmetry due to finite N s is the minimal 
provided that the weights {uj s } are fixed according to the formula derived in 
Ref. [6]. Note that in this framework, the bare mass m q (no matter heavy 
or light) in the valence quark propagator (D c + m g ) _1 is well-defined for any 
gauge configurations. 

We generate 100 gauge configurations with single plaquette gauge action 
at (3 = 7.2 on the 32 3 x 60 lattice. For m = 1.3 and N s = 128, we fix the 
weights {u s } with X min = 0.1 and X max = 6.4, where \ min < \{\H W \) < X max 
for all gauge configurations. For each configuration, point to point valence 
quark propagators are computed for 33 bare quark masses in the range 0.01 < 
m q a < 0.85, with stopping criteria 10" 11 for the conjugate gradient. Then the 
norm of the residual vector of each column of the quark propagator is less than 
2 x 10" 11 



\\(D c + m q )Y - I|| < 2 x 10" 
and the chiral symmetry breaking due to finite N s is less than 10~ 14 , 



a 



ytQ2 y 
1 °opt 1 



< lO" 14 , 



In this paper, we measure the time-correlation function 

C T (t) = (^tr{T(D c + m Q )7 x } i T(D c + m q )^ (2) 

for scalar (S), pseudoscalar (P), vector (V), axial- vector (A), and tensor (T) 
mesons, with Dirac matrix T = {1, 75, 7i, 757*, 757474} respectively. For vector 
mesons, we average over i — 1, 2, 3 components, namely, 

Similarly, we perform the same averaging for axial-vector and tensor mesons. 

The time-correlation function Cr{t) is measured for the following three 
categories: 

• Symmetric masses with tuq = m q for 33 quark masses. 

• Asymmetric masses with fixed mq = = 0.68a -1 , and m q running over 
all different quark masses. 

• Asymmetric masses with fixed tuq = m c = 0.16a _1 , and m q running over 
all different quark masses. 
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3 Determination of a \ nib, win an d m s 

In Ref. [8], we determine the inverse lattice spacing from the pion decay 
constant, with experimental input f n = 131 MeV. However, in this paper, we 
do not use the same method since the smallest quark mass turns out to be 
rather heavy (~ m s /2), thus chiral extrapolation to m q ~ does not seem 
to be feasible. Nevertheless, we can use the mass and decay constant of the 
pseudoscalar meson ?7 C (2980) to determine m c and a -1 simultaneously. This 
can be seen as follows. 

For symmetric masses m,Q = m q , the pseudoscalar time-correlation function 
Cp{t) (r = 75) is measured, and fitted to the usual formula 

\c~ mpat + e~ mpa ( T ~ t h 

2m pa 

to extract the mass rripa and the decay constant 

f P a = 2m q a I . 

mpa^ 

Then the ratio rrip/fp can be obtained for each m q . 

Since f Vc has not been measured in HEP experiments, we do not have 
a physical value for the ratio m Vc /f Vc . Nevertheless, we can obtain the ra- 
tio m Vc / j nc = 6.8(2) from our previous study of pseudoscalar mesons on the 
20 3 x 40 lattice at /3 = 6.1 [8]. Then we can use this theoretical value 6.8 
to discriminate which m q can give the ratio mp/fp closest to this value. We 
find that at m q a = 0.16 the ratio mp/fp = 6.8(1), which is the closest to 6.8. 
Thus we fix m c a = 0.16. Then we use the experimental mass of ?7 C (2980) to 
determine a -1 through the relation 

m P a\ mc = (2980 MeV) x a = 0.388(3) 

and obtain a -1 = 7680(59) MeV. To check the goodness of the values of m c and 
a" 1 , we compute the time-correlation function of cr/jC, and extract the mass of 
the vector meson equal to 3091(11) MeV, in good agreement with J/\I r (3097). 

The bare mass of strange quark is determined by extracting the mass of 
vector meson from the time-correlation function Cy(£). At m q a = 0.02, mya = 
0.1337(5), which gives my = 1027(38) MeV, in good agreement with the mass 
of ^>(1020). Thus we take the strange quark bare mass to be m s a = 0.02. 
Similarly, at m q a = 0.68, mya = 1.2308(4), which gives my = 9453(3) MeV, 
in good agreement with the mass of T(9460). Thus, we fix the bottom quark 
bare mass to be m^a = 0.68. 

Note that the spatial size of our lattice (L ~ 0.8 fm) seems to be small 
at first glance, however, even for the smallest quark mass m q a = 0.01, its 
pseudoscalar mass satisfies mpL > 4, thus the finite size effects should be well 
under control. 
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4 Charmonium cc and Bottomonium bb 

First of all, we check to what extent we can reproduce the charmonium mass 
spectra which have been measured precisely by high energy experiments. 

Our results of the mass spectra of the lowest-lying states of charmonium 
are summarized in Table 1. The first column is the Dirac matrix used for 
computing the time-correlation function (2). The second column is J PC of 
the state. The third column is the conventional spectroscopic notation. The 
fourth column is the \t m i n ,t max ] used for fitting the data of Cp(t) to the usual 
formula 

\ P ~Mat , -Ma(T-t)] 



to extract the meson mass M and the decay amplitude z. The fifth column is 
the mass M of the state, where the first error is statistical, and the second is 
our estimate of systematic error based on all fittings satisfying x 2 /dof < 1.3 
and \t max — t min \ > 6 with t min > 10 and t max < 50. The last column is the 
corresponding state in high energy experiments, with the PDG mass value [9]. 
Evidently, our mass spectra of charmonium are in good agreement with the 
PDG values. Note that our result of the hyperfine splitting (l 3 5i — l 1 ^) is 
111(14)(18) MeV, comparing with the PDG value 118 MeV. 

For the pseudoscalar r] c , we also obtain its decay constant f Vc together with 
its mass, through the decay amplitude z in the equation 



f„ c a = 2m c a 

mi a z 



Our result is 



f Vc = 438 ± 5 ± 6 MeV (3) 

where the first error is statistical, and the second is our estimate of systematic 
error based on all fittings satisfying x 2 /dof < 1.3 and \t max — t min \ > 6 with 
tmin > 10 and t max < 50. So far, f Vc has not been determined in high energy 
experiments. 

Next, we turn to the bottomonium (bb) states. Our results of the mass 
spectra of the lowest-lying states of bottomonium are summarized in Table 2. 

First, we look at the pseudoscalar rj b . It was first reported by ALEPH 
Collaboration [13]. However, it has not been confirmed by other HEP exper- 
imental groups. Thus we suspect that its mass 9300(20) (20) MeV might not 
have been determined accurately. It is interesting to see whether the mass of 
rjb will turn out to agree with our theoretical value 9380 ±10 MeV. Besides 
the mass of rjb, we also determined its decay constant 

f Vb = 801 ± 7 ± 5 MeV (4) 
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Table 1: The mass spectra of lowest-lying (n = 1) charmonium cTc states 
obtained in this work, in comparison with the PDG values in the last column. 



r 


jPC 


n 2S+1 Lj 




X 2 /dof 


Mass ( Me V) 


PDG 


i 


0++ 


l 3 Po 


[19,42] 


0.32 


3413(14)(9) 


Xco(3415) 


75 


o-+ 


I 1 So 


[22,38] 


1.02 


2980(10)(12) 


r/ c (2980) 


li 


\— 


1 3 S 1 


[19,38] 


1.18 


3091(11)(14) 


J/^(3097) 


7s7i 


1++ 


l 3 Pi 


[18,42] 


0.51 


3516(13)(8) 


Xci(3510) 


75747i 


1+- 




[19,43] 


0.41 


3526(13)(9) 


/i c (3524) 



Table 2: The mass spectra of lowest-lying (n = 1) bottomonium brb states 
obtained in this work. The last column is the experimental state we have 
identified, and its PDG mass value. 



r 


jPC 


n 2S+1 Lj 




X 2 /dof 


Mass(MeV) 


PDG 


i 


0++ 


l 3 P 


[21,39] 


0.21 


9863(15)(8) 


X6o(9859) 


75 


o- + 


I 1 S 


[27,35] 


0.72 


9383(4) (2) 


^(9300) ? 


li 


1" 


l 3 5i 


[20,39] 


1.19 


9453(3) (2) 


T(9460) 


75% 


1++ 


1 3 P! 


[22,38] 


0.13 


9896(20) (8) 


X6i(9893) 


75747* 


1+- 




[22,38] 


0.10 


9916(22)(8) 





where the first error is statistical, and the second is our estimate of systematic 
error based on all fittings satisfying x 2 /dof < 1.3 and \t max — t min \ > 6 with 
tmin > 10 and t 

in ax 

< 50. 

Finally, we note that the tensor meson hb has not been observed in high 
energy experiments, thus our result of its mass 9916 ± 30 MeV serves as the 
first prediction from lattice QCD with exact chiral symmetry. Even though it 
is obtained in the quenched approximation, we suspect that it might provide 
a reliable prediction of hadron mass spectra, especially for mesons with heavy 
b quark. 

5 Mesons with quark contents sb and cb 

The decay constants of heavy- light pseudoscalar mesons (e.g., fs s , /d, and 
fo s ) play an important role in extracting the CKM matrix elements which are 
crucial for testing the Standard Model via the unitarity of CKM matrix. The- 
oretically, lattice QCD with exact chiral symmetry provides a reliable frame- 
work to compute the masses and decay constants of heavy-light pseudoscalar 
mesons nonperturbatively from the first principles of QCD. Our theoretical 
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predictions of fo and fo a have been presented in Ref. [8], which turn out 
in good agreement with the recent experimental results [10, 11] from CLEO 
Collaboration. For a recent review of heavy-light pseudoscalar decay constants 
from lattice QCD, see Ref. [12] and references therein. 

Before we present our results of f Bs and m^, we recall the basic formulas 
as follows. In general, the decay constant fp of a heavy-light pseudoscalar 
meson P is defined as 

(0\A„(0)\P(f)) = ip„f P 
where A^x) = Q{x)^f tl r )^q{x) is the axial-vector current. Using the formula 

d^A^ = (m q + m Q )Qj 5 q 

one obtains 

|<o|g T 5s|P(o)>| 



fp = (m q + m Q )- 



m. 



where the pseudoscalar mass m P a and the decay amplitude z = \ (0\Q^5q\P(0))\ 
can be obtained by fitting the pseudoscalar time-correlation function Cp(t) to 
the usual formula 



7 2 

\ e -m P at _|_ e -m P a(T-t)-\ 



2mpa 

Our results for B s are: 

m Bs = 5385 ± 27 ± 17 MeV, f Bs = 253 ± 8 ± 7 MeV (5) 

where the first error is statistical, and the second is our estimate of systematic 
error based on all fittings satisfying x 2 /dof < 1.3 and \t max — t min \ > 6 with 
tmin > 10 and t max < 50. Our result of mp a is in good agreement with 
the experimental value 5368 MeV compiled by PDG. Since fp a has not been 
measured in high energy experiments, our result serves as the first prediction 
from lattice QCD with exact chiral symmetry. 

Next, we present our results of the mass spectra of the lowest-lying states 
of beauty mesons with quark content sb, which are summarized in Table 3. 
Here we have identified the scalar bs meson with the state _B*j(5850) observed 
in HEP experiments, due to the proximity of their masses. Theoretically, 
this implies that -B*j (5850) possesses J p = + , which can be verified by HEP 
experiments in the future. Moreover, we have obtained the masses of the axial- 
vector and tensor mesons which have not been observed in HEP experiments, 
clS db prediction from lattice QCD. 

Next, we turn to the heavy mesons with beauty and charm. For the pseu- 
doscalar B c , we obtain 

m Bc = 6278 ± 6 ± 4 MeV, f Bc = 489 ± 4 ± 3 MeV (6) 
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Table 3: The mass spectra of lowest-lying brs meson states obtained in this 
work. The last column is the experimental state we have identified, and its 
PDG mass value. 



r 


J p 


n 25+1 Lj 


\tmim tmax\ 


X 2 /dof 


Mass ( Me V) 


PDG 


i 


0+ 


l 3 Po 


[20,40] 


0.38 


5852(15)(12) 


B*j(5850) 


75 


0" 


1% 


[26,32] 


0.56 


5385(27)(17) 


P s (5368) 


7i 


1- 


1 3 S 1 


[25,33] 


0.50 


5424(28)(19) 


B*(M12) 


757i 


1+ 


l 3 Px 


[18,42] 


0.68 


5884(16)(13) 




75747i 


1+ 




[18,38] 


0.62 


5897(16)(12) 





Table 4: The mass spectra of lowest-lying brc meson states obtained in this 
work. The last column is the experimental state we have identified, and its 
PDG mass value. 



r 


J p 


n 2S+1 Lj 




X 2 /dof 


Mass(MeV) 


PDG 


i 


0+ 


1 3 P 


[19,41] 


0.26 


6732(13)(9) 




75 


o- 


1% 


[19,38] 


1.14 


6278(6) (4) 


P c (6286) 


7i 


1- 


1 3 S 1 


[19,38] 


1.31 


6315(6)(5) 




757i 


1+ 


l 3 Pi 


[19,41] 


0.27 


6778(12)(7) 




75747i 


1+ 




[18,42] 


0.45 


6796(10)(7) 





where the first error is statistical, and the second is our estimate of systematic 
error based on all fittings satisfying x 2 /dof < 1.3 and \t max — t min \ > 6 with 
tmin > 10 and t max < 50. Our result of m Bc is in good agreement with 
the experimental value 6286(5) MeV measured by CDF Collaboration [14], 
while fs c has not been measured by HEP experiements. In principle, fs c can 
be measured from the leptonic decay — > since its decay width is 

proportional to f^jVcbl 2 - 

In Table 4, we summarize our results of the mass spectra of the lowest- 
lying states of mesons with beauty and charm. Except for the pseudoscalar 
B c , other states have not been observed in experiments. It will be interesting 
to see to what extent the experimental results would agree with our theoretical 
values. 

6 Masses of b and c quarks 

In the Standard Model, the quark masses are fundamental parameters which 
have to be determined from high energy experiments. However, they cannot 
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Table 5: The decay constants of pseudoscalar mesons obtained in this work, 
together with their masses. They are identified with the corresponding PDG 
mesons listed in the last column, however, the decay constants have not been 
measured in HEP experiments. 



Qr<? 




xYdof 


Mass(MeV) 


/p(MeV) 


PDG 


b7 5 b 


[27,35] 


0.72 


9383(4)(2) 


801(7)(5) 


r/ fe (9300) 


b7 5 c 


[19,38] 


1.14 


6278(6) (4) 


489(4) (3) 


5 C (6287) 


b7 5 s 


[26,32] 


0.56 


5385(27)(17) 


253(8)(7) 


5,(5368) 


C75C 


[22,38] 


1.02 


2980(10)(12) 


438(5)(6) 


7^(2980) 



be measured directly since quarks are confined inside hardrons, unlike an iso- 
lated electron whose mass and charge both can be measured directly from its 
responses in electric and magnetic fields. Therefore, the quark masses can only 
be determined by comparing a theoretical calculation of physical observables 
(e.g., hadron masses) with the experimental values. Evidently, for any field 
theoretic calculation, the quark masses depend on the regularization, as well 
as the renormalization scheme and scale. One of the objectives of lattice QCD 
is to compute the hadron masses nonperturbatively from the first principles, 
and from which the quark masses are determined. 

We have used the mass of the vector meson T(9460) to fix the bare mass 
of b quark equal to m b = 0.68a -1 = 5.22(4) GeV 1 To transcribe the bare 
mass to the corresponding value in the usual renormalization scheme MS in 
high energy phenomenology, one needs to compute the lattice renormalization 
constant Z m = Z~ x , where Z s is the renormalization constant for ipip. In 
general, Z m should be determined nonperturbatively. However, in this case, 
the lattice spacing is rather small (a ~ 0.026 fm), thus we suspect that the 
one-loop perturbation formula [15] 



a 2 r 1 
ZAfj) = 1 + f- ln(aV) + 0.17154 , (m = 1.30) (7) 
47H L J 

already provides a very good approximation for Z s . At (3 = 7.2, a -1 = 
7.680(59) GeV, and fi = 2 GeV, (7) gives Z s = 1.086(2), which transcribes 
the bare mass m b = 5.22(4) GeV to 



m^ S (2 GeV) = 4.81 ± 0.04 GeV 



where the error is due to the uncertainty in the lattice spacing. Now if we also 
include the uncertainty in the determination of the b quark bare mass, which 



1 Note that we have not tuned the value of nib, thus the mass of T is 9453(3)(2) MeV 
(see Table 2) rather than 9460 MeV. 
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is estimated to be 5{m b a) = 0.002, then we obtain 



m^ S (2 GeV) = 4.81 ± 0.05 GeV (8) 
In order to compare our result with the PDG average, we have to obtain 



at the scale /i = m;,. This can be obtained by solving the equation 
rh b = m b Z m (fi = m b ). Our result is 



mj^ S (m 6 ) = 4.65 ± 0.05 GeV (9) 

which seems to be higher than the PDG average 4.20 ± 0.07 [9]. 

Now we turn to the c quark mass. Using (7), the c quark bare mass 
m c = 0.160(5)a~ 1 is transcribed to 



m^ S (2 GeV) = 1.13 ± 0.04 GeV (10) 

where the error incorporates the uncertainties in the lattice spacing and the 
c quark bare mass. Again, to compare our result with the PDG average, we 

obtain at the scale /i = m c , 



m^ S (m c ) = 1.16 ±0.04 GeV (IT 
which is in good agreement with the PDG average 1.25 ± 0.09 [9]. 



7 Concluding remark 

We have performed the first study of treating b, c, and s quarks as Dirac 
fermions in lattice QCD with exact chiral symmetry. The lowest-lying mass 
spectra of mesons with quark contents bb, cb, sb, and cc are determined, 
together with the pseudoscalar decay constants. Our results of the meson 
mass spectra in Tables 1-4, pseudoscalar decay constants in Table 5, and the 
b and c quark masses in (9) and (11), suggest that lattice QCD with exact 
chiral symmetry is a viable framework to study heavy quark physics from the 
first principles of QCD. 

For systems involving u/d quarks, one may use several quark masses in the 
range < m q < m s to perform the chiral extrapolation. To this end, one 
may choose a coarser lattice (e.g. (3 = 7.0), then it is possible to accommodate 
a wide range of quark masses m s /A < m q < m b on the 42 3 x 64 lattice, without 
significant discretization and finite-size errors. This study is now in progress. 
Evidently, it has become feasible to treat heavy and light quarks as Dirac 
fermions on the lattice, in lattice QCD with exact chiral symmetry. 
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